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Abstract 

We show that the Polyakov loop approaches its asymptotic value at high temperatures from 
below in any strongly-coupled gauge theory at large N c with a dual supergravity description. In 
these theories the Polyakov loop can be used to define a quantity that behaves as a free energy of 
a heavy quark in the medium. This should be compared to the recent lattice results in pure glue 
SU(3) theory at high temperatures where the Polyakov loop approaches its non-interacting limit 
from above. We show that the "heavy quark free energy" obtained from the renormalized loop 
computed on the lattice does not behave as a true thermodynamic free energy. We conjecture that 
this should always be the case in asymptotically free gauge theories. 
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SU(N C ) gauge theories considerably simplify when the number of colors, N c , is large 
The deconfmement phase transition in this limit is predicted to be a strong first order 
transition where the latent heat ~ iV| and, in fact, this has been shown to be true in 
lattice calculations performed in pure gauge theories with N c > 3 |2J. In this paper I show 
that when N c — > oo the expectation value of the (regularized) Polyakov loop, which is the 
order parameter [3| for the deconfining phase transition in SU (N c ) gauge theories without 
dynamical fermions, becomes a monotonic function of the temperature in the deconfined 
phase of gauge theories with supergravity duals. 

The order parameter for deconfinement in an SU(N <i ) theory without dynamical fermions 



is defined in terms of the path ordered Polyakov loop [3] 

L(f) = Pe ! ^o 1/Tio(S ' T)dT (1) 

where A is the non-Abelian gauge field operator and g is the coupling. This operator 
becomes gauge invariant after performing the trace over the fundamental representation of 
SU(N C ), denoted here as tr. Note that trL transforms as a field of charge one under the 
global Z(N C ) symmetry, i.e., trL — > e i27ra//iVc tr L where a = 0, . . . , N c — 1. The magnitude 
of the expectation value of this operator after the (normalized) trace is denoted here as 
£ = |(tr L)|/iV c . Below T c the system is Z(N C ) symmetric, which implies that (trL) = 0. 
Above T c this global symmetry is spontaneously broken, £(T) ^ 0, and the system lands in 
one of the possible Z(N C ) vacua. 



Wilson lines are renormalizable operators [4] . The expectation value of a Wilson line has 
a mass divergence that depends on the length of the loop |4j and for a thermal Wilson line the 
length is 1/T and £ ~ exp(— m dlv /T), where m dlv is a divergent mass that is linear in the UV 
cutoff. On a lattice with spacing a one would have m dtv ~ 1/a. It is expected that m dlv > 
and, thus, the bare quantity computed on the lattice vanishes in the continuum limit jsj. 
The wave function renormalization of the loop is given in terms of Z = exp(—m dm /T), which 
defines then the renormalized quantity L^ = L/Z [5( (we only consider smooth loops here). 
We shall assume that all the quantities discussed here have been properly renormalized so 
that we can drop the subscript R in the following. Moreover, note that £ is a positive-definite 
number that vanishes below T c in pure glue gauge theories. 

The thermal average of the 2-point function C(r,T) = (trLt(r) trL(0))/A r ^ is usually 
associated with the free energy of an infinitely heavy QQ pair separated by a distance r 



2 



in the medium 



It is appropriate at this point to re-examine this argument when N c is 



large. We use the standard definition for the thermal expectation value of a given operator 



B 
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Tr B e~ H/T 



where % is the Hamiltonian of the system and Tr denotes the trace operation defined in 
terms of a path integral over the fields. 

The free energy of a system composed by gluons and a QQ pair is 

F(n, f 2 , T) = — T In J2{s\e- H/T \s) (3) 

\s) 

where TL is the Hamiltonian of the whole system (including the quarks), and the sum is over 
all the states \s) with a heavy quark at r\ and an anti- heavy quark at r 2 (with r = |r*i — r 2 |). 
After introducing creation and annihilation operators for the quarks and anti-quarks and 
regularizing the loops, one can define a renormalized QQ interaction energy 

AF QQ (r,T) = F(r,T)-F\T) = -T\nC(r,T) (4) 

in terms of the difference between the free energy of the system with the quarks, F(r,T), 
and the free energy of pure glue at finite T, F 9 {T). 

The total free energy of the pure glue medium can be written as 

F9{T) = N 2 C F°{T) + F 9 (T) + 0{1/N 2 C ) (5) 

where the coefficients Ff = F[(X,T) are only functions of T and the t' Hooft coupling 
AyM = 9ymN c - A similar expansion must also hold for the system which includes the effects 
from heavy quark probes 

F(r, T) = N 2 C F 2 (T) + F (r, T) + 0(1/N 2 C ) (6) 

where again Fi = Fi(X,T) but, in general, Fi and Ff are not identical. However, note that 
F2 — F2 since at large N c the contribution from the quarks cannot enter at that order. 
Therefore, one can see that 

hm AFqq (r, T) = F (r, T) - F 9 (T) . (7) 

iv c — >-oo 

All the coefficients in the series in Eqs. (jSJ) and ([6]) define a well-defined free energy order 
by order in N c , i.e, Fi is the contribution of order N l to the total free energy, which must 



behave as true thermodynamic free energy 9J as well. This argument is certainly valid when 
N c — > oo and all the coefficients in the series above can be unambiguously defined. 

Note that £ 2 (T) = C(r — > oo, T) and with this relation one can obtain a quantity that is 
associated with a single heavy quark in the medium via 

F Q {T) = -Thx£{T). (8) 

The standard physical interpretation is that difference of free energies AF^^r, T) repre- 
sents the free energy contribution of the heavy QQ pair to the medium 6|. In general, the 
quantity defined in Eq. (J3J) may not behave like a free energy in the thermodynamic sense. 
The difference between free energies is not guaranteed to be a free energy since the static 
color-electric field sourced by the QQ pair can considerably change the energy of the system. 
For instance, in a confining theory the energy of the system below T c is not a well defined 
quantity if the quarks are infinitely distant from each other. 

Let us now see under what conditions Fq(T) (and hence AFgg(r, T)) can be associated 
with a free energy. The derivative of I with respect to T is 

where Uo(r, T) = Fq(t, T) — TdFo(r, T) / dT is the internal energy of the system that contains 
a quark while Uq{T) = Fq{T) — T dF$ (T) / dT is the corresponding quantity obtained in the 
absence of a quark. As argued above, each of these internal energies should display the 
expected thermodynamic behavior when defined in terms of the free energy 9]. 

Since in a confining theory I is zero at T < T c and positive otherwise, we see that 
immediately above T c the derivative dtjdT > 0. This means that the difference in internal 
energies in Eq. ([9]) is positive at T — > T c + . In homogenous systems in equilibrium the internal 
energy can only increase with T since the specific heat must be non-negative 9j. However, 
this does not imply that di/dT should be always positive. In fact, one can imagine that 
because of the quark the specific heat of the system in this case can increase at a different 
rate and at higher temperatures the difference in internal energies in the equation above can 
vanish and then become negative. If this happens, Fq(T) cannot be a free energy because 
its corresponding "internal energy" Uq(T) = Uo(r — > oo,T) — Uq(T) would have negative 
specific heat Cq(T) = dUq(T) / dT < 0. On the other hand, if the derivative of the Polyakov 
loop remains positive at any finite T > T c then Fq(T) acquires a true thermodynamic 



meaning, Cq(T) > 0, and d£(T)/dT — > + when T — > oo where there is usually an UV fixed 
point. 

Perturbative QCD calculations in the hard-thermal-loop approximation performed using 



a running coupling constant 



cently, lattice simulations 
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11] suggest that I should approach one from above. Re- 



13] in SU (3) pure glue showed that t increases very rapidly 



after its jump at T c and, in fact, it overshoots one around 3T C and keeps increasing until 
it reaches a maximum near 10T C [kJ where it starts to decrease with increasing temper- 
atures towards a constant value that is close to one. According to the discussion above, 
the fact that in these calculations d£/dT changes sign shows that Fq(T) (or equivalently, 
AFqq(t, T)) cannot be identified with a quark free energy in the QCD plasma. In this case, 
the associated "heavy quark specific heat" cq < 0. The thermal Wilson line is just the trace 
of the propagator of a quark which, because of its infinite mass, moves in a straight line in 
imaginary time [14| . This quantity obviously depends on T but this dependence need not 



be that of a thermodynamic function such as the free energy. 

This issue regarding the physical interpretation of the Polyakov loop correlator was dis- 
cussed, for instance, in jljj while in [ijj] it was shown that in the presence of a finite baryonic 
chemical potential Fq is expected to increase with /i, which clearly shows that this quantity 
is not a true free energy. It is more appropriate perhaps to call Fq(T) simply the heavy 
quark interaction energy. 

The dominant finite T contribution to the binding energies of very tightly bound pairs, 
such as the Y , comes from Fq(T) 16Ml8|. If in QCD this quantity is not a true free energy, 
then the usual "entropy subtraction" performed in quarkonia binding energy studies 16] in 
pure glue theories may be less justified. This will affect quantitative analyses of quarkonia 
melting, which is an important signature for deconfinement (if] ]. 

It is natural to think that asymptotic freedom may be the reason for this non-monotonic 
behavior with T displayed by the Polyakov loop. In fact, the hard-thermal-loop calculation 



of I performed at leading order gives I ~ 1 + Ayjvr(m£)/T), where mo ~ V^ym T is the 
Debye mass. For a fixed value of the coupling m^/T is a constant and i is a number larger 
than one. Clearly, once the coupling is allowed to run dijdT < because dXyu/dT < 0. 
Thus, while confinement implies that immediately above T c the derivative of i is positive, 
at high T Debye screening and asymptotic freedom impose that dijdT must be negative. 
These general arguments lead us to conjecture that in confining gauge theories with a trivial 
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UV fixed point there must be a finite value of temperature, say T*, where d£(T*)/dT = 0. 
The exact value of T* may not be obtained perturbatively unless it is sufficiently larger than 
T c . As discussed above, in this class of theories Fq(T) cannot be a true thermodynamic free 
energy. In fact, since d£/dT < at high T this implies that the specific heat of the gluon 
medium decreases (by a small amount) due to the presence of the heavy quark. 

Since the thermodynamic properties of pure gauge theories computed on the lattice do 
not change significantly when N c > 3 Q], one should expect that £ is not always a monotonic 
function of T in QCD at large N c . The published lattice data for the Polyakok loop in full 
SU (3) QCD with dynamic fermions fjj] is still restricted to a much narrower temperature 
range than the calculations performed in pure glue 13] . The general arguments discussed 
here indicate that in full QCD £ is not be a monotonic function of T, which means that 
Fq should not be associated with a free energy. It would be interesting to see if T* in this 
case shifts towards the transition region or remains close to the value in pure glue ~ 10T C 
[l^j ]. However, one should also expect that the exact location of T* could depend on the 
regularization scheme used to obtain to compute £ on the lattice. 

What happens in gauge theories with nontrivial UV fixed points? Since these theories 
are not asymptotically free, one may expect that the general argument presented above does 
not apply and in this case Fq(T) may behave as a free energy. The gauge/string duality 



3. 



2l| provides a way to study the behavior of the Polyakov loop in this class of theories. 



We shall show in the following that in confining gauge theories with a supergravity-like dual 
the Polyakov loop always approaches its asymptotic value at high T from below. 

Let us first recall how the deconfinement phase transition takes place in these theories 



2l| . The general idea is that the deconfined phase of the 4d gauge theory is dual to 



a theory of gravity defined in a higher-dimensional geometry with a black brane whose 
Hawking temperature equals the temperature of the plasma. The free energy of the system 
above T c equals the free energy of the black brane and is of order AT? (because the effective 
gravity coupling ~ while at lower temperatures another solution of the supergravity 

equations that does not have a horizon has a smaller free energy, which is of order iV° after 
the inclusion of quantum corrections. The deconfinement phase transition in this approach 
is predicted to be a strong first order transition with latent heat ~ N%, which agrees with 
lattice calculations of pure glue gauge theories at large N c [2]. 

We shall assume in the following that N c — > 00 and the radius of the background space- 
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time R is much larger than the string length \fa' (the supergravity limit). For instance, in 
the planar limit of N = 4 Supersymmetric Yang-Mills (SYM) theory the t'Hooft coupling 
is Xsyai — R A /<y' 2 20|, |2l|. When a theory admits a supergravity description several sim- 
plifications occur. For example, because the classical supergravity action goes like ~ N%, 
the Euclidean path integral is dominated by the spacetime that is a solution of the classical 
gravitational equations of motion and it becomes ~ exp(— S sugra ) with the on-shell action 



being S sugra ~ N% J d w x \/—G (1Z + . . .) while . . . denote the contribution from the other 
supergravity modes. 

An infinitely massive quark in the 4d non-Abelian gauge theory at finite T is dual to 
a string in the bulk that hangs down from a probe D-brane at infinity and reaches the 
black hole horizon in the bulk spacetime 22j . The holographic representation of £ is simply 
£{T) = exp(— 5jv<?(^)) where Sng is the classical Nambu-Goto action (in Euclidean time), 
which is proportional to R 2 /a', and D is the worldsheet that is a solution of the equations 
of motion describing a single string in the bulk (see 23|-|25|). Here the string is indeed a 
real probe that does not backreact on the geometry and the full partition function factorizes 
Zq(T)/Z(T) = exp(— Sng) (we are assuming that the divergence in the string action that 
appears near the boundary has been already properly removed). 

Similar arguments to those used in Eqs. (jSHTJ) also hold in this case but the important 
difference here is that F 9 (T) = N^F^T) when N c — > oo and R 2 /a' — > oo, i.e, the other 
coefficients in the N c expansion of the theory without the quark probe vanish in this clas- 
sical approximation. Also, should not depend on R 2 /a' in this limit 26]. This means 
that in the classical supergravity approximation Fq(T) = Fo(r — > oo,T)/2. Thus, in the 
supergravity limit Fq(T) behaves as a free energy and, according to the general discussion 
above, I should be a monotonic function of T. The free energy in this limit is simply given 
by the Nambu-Goto action via Fq = T Sng evaluated using the background spacetime. 
The specific heat of the plasma increases (by a very small amount) when the heavy quark is 
added, which is the opposite of what occurs in a pure glue plasma according to the discussion 
presented earlier involving the lattice data. 



Since screening also takes place at strong coupling [23], it seems that the reason why £ 
behaves this way in these theories is because they are not asymptotically free. One may 
argue that holographic methods can provide an approximation for £ in QCD only when 
T < T* since above T* effects from asymptotic freedom are not negligible (even though the 
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coupling around T* may not be small). However, when corrections to supergravity of order 
R 2 / (a' N 2 ) are present in the bulk the string cannot be seen as probe and it may be possible 
to show that there is a finite T* where di/dT vanishes. 

Thus, we see that the temperature dependence of I in confining theories at large N c 
varies with the type of fixed point present in the UV. For theories that are asymptotically 
free there must be a finite value of T > T c where d£(T)/dT = while for theories with 
nontrivial fixed points described by supergravity duals i is a monotonic function of T in 
the deconfined phase. However, we would like to stress again that once corrections to the 
supergravity background are taken into account i may not be a monotonic function of T. 
This could be the case even if this higher order theory is not asymptotically free. 

This result for the T dependence of £ in gauge theories computed in the supergravity 
approximation can be used to show that strongly-coupled gauge theories at large N c with 
nontrivial UV fixed points may display unphysical results depending on the scaling dimen- 
sion, A < 4, of the operator used to deform the original CFT and induce a nontrivial 
renormalization group flow. The general class of gravity duals we are interested in were 
studied in detail in [27| where it was shown that it is possible to obtain black brane back- 
grounds with equations of state that are similar to the QCD equation of state obtained on 
the lattice (see also [28] for a model involving the dilaton). 

The gravity action describes the interactions between a 5d metric (in the Einstein frame) 



and a single scalar field (f> 27] 



A ~ 7^2 I '/'"'-''V -('< 



2 



*-^-v<#) 



(10) 



'5 

where k\ ~ 1/N 2 . Conformal invariance in the UV can be obtained by imposing that — > 
at the boundary and 

1 2 1 

lim V(0) = ~ + ^A(A - 4)0 2 + O(0 4 ), (11) 
where the mass squared of the scalar m 2 ,R 2 = A(A — 4) > —4 and, according to the 



Breitenlohner-Freedman bound [29], 1 < A < 4. Eq. ( TTTT) implies that the spacetime will be 
asymptotically AdSs with radius R. The deviations from conformal behavior only become 
significant deep down in the geometry when u is small. The boundary of the bulk spacetime 
is located at u — > oo. The Ansatze for the Einsten frame metric and the scalar are 

dv 2 

ds 2 = a 2 (u)(-f(u)dt 2 + dx 2 ) + p,™ f ,y </> = 0(«) • (12) 



Finite temperature effects are included by considering solutions that display a horizon at 
u = Uh where f(uh) = 0. At temperatures T/A ^> 1, where A is the energy scale associated 
with the scalar, one should be close to the AdS 5 solution and in this case (fih = <f)(uh) <C 1. 
A consistent high T expansion of several transport coefficients in gravity duals with the 



action (fTU|) and 2 < A < 4 was developed in 30] and a generalization of these results when 



1 < A < 4 was done in [31(. At temperatures where (fih *C 1 the corrections to the CFT 
values (the T — > oo limit) start at <p 2 h where 

< r(A/4)! m 



(A - 2lr(A/2 — 1) \nT J 



25j, one can see that Nambu-Goto action or Fq(T) is 



Following the results presented in 
given by 

F Q (T)=T S £ duq^u))^ (14) 

apart from a temperature independent term used to remove the divergence at u — > oo. Here 
T s = l/(27ra') is the fundamental string tension and q(<fi) is a function that describes how 
the fundamental string couples to the scalar. In the supergravity approximation to type IIB 
string theory, assuming that is the dilaton this coupling function in the Nambu-Goto action 
would be equal to exp(0/2). M = 4 SYM is conformal invariant and in the supergravity 
approximation the dilaton is a constant. In this case the gauge coupling Xsym does not 



run and the expression for the free energy is Fq = —yXsYM T/2 at strong coupling (the 
linear dependence with T follows from conformal invariance). Note that in exactly conformal 
theories such as M = 4 SYM the internal energy always vanishes, i.e, Uq = to all orders 
in the t'Hooft coupling. 

In the case of the effective model in Eq. ( 1101) the function q(<ft) can be in principle an 
arbitrary function of 0. However, we shall show in the following that it is possible to obtain 
constraints on A and q(<ft) by imposing that Fq in Eq. ffT4l) behaves as a free energy (or, 
equivalently, that £ is a monotonic function of T) in gauge theories dual to supergravity 
backgrounds that have a nontrivial fixed point in the UV. These constraints are derived 
using the high T expansion of Uq and Cq computed from Eq. (fT4l) . Since in these gravity 
duals 

and Uq(T ~ T c ) > 0, the condition d£/dT > implies that Uq should remain positive at 
any T and, thus, Cq > 0. While near the transition these quantities can only be computed 



numerically, one can use the expansion defined in 30j to determine their high T behavior 
analytically. In fact, A/T < 1 implies that <ph *C 1 and one can expand the coupling 
function in a power series around (ph ~ 0, q(4>h) = ?(0) + q'(0)(ph + ^"(0)0^/2 + . . ., to show 
that the heavy quark entropy is 

n=0 ; 



where 



and the first coefficients are 



rf < A > - (A - 2)^/2 -1) (17) 



and 



So (A) = l, Sl (A) = ^, (18) 

«a(A) = |^|| + ^(4 - A)(2 - A) tan(vrA/4) + ■ (19) 

It is easy to show using Eq. ( JT6l) that when A^3 the internal energy is 

p2 ( A — A")n / A \ (4-A)n 

W - ^(°) r E «.(A)i(A)- (ir) ' (2 °) 

which leads to the following expression for the specific heat 

d2 /AN ( 4 ~ A ) n 

C (r) = ^(O) £ *„(A)d(A)» (A - 4)n — . (21) 

n=l N^ ' 

Note that Eq. fl20l) implies that dljdT — >■ when T/A — > oo. Let us assume that the coupling 
function is such that q(0) > and si(A) 7^ 0. Then, by requiring that Cq is positive at 
high T one obtains that si(A) < (note that d(A) > when 1 < A < 4). However, if 
si(A) is negative then Uq can only be positive if A > 3 (the generalization of the formulas 
above for the case where A = 3 is straightforward). One can see from Eq. ( TT8|) that this 
condition implies that q'(0) < 0. Thus, simple coupling functions that increase exponentially 
with lead to inconsistencies in the supergravity approximation of gravity theories dual to 
CFTs deformed by a relevant operator (if A = 4 the standard exponential coupling is, of 
course, valid). Exponentially decreasing coupling functions, such as exp(— with £ > 0, 
usually appear in the description of Dl-branes whose endpoints at the boundary correspond 
to magnetic monopoles in the plasma. Thus, since the condition q'(0) < cannot probably 
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be easily understood in terms of fundamental strings, one may want to consider the case 
where q'(0) = 0. If si(A) is identically zero then 52(A) must be negative so that Cq > 0. 
This implies that Uq is positive at large T only if A > 3.5. Note also that this constraints 
the value of q"(0) in terms of A via Eq. (TTTJT) . Perhaps a good choice for phenomenological 
applications would be q(<p) = cosh(£0), where £ = £(A) is chosen to give s 2 (A) < 0. This 
coupling function becomes exponentially large in the IR and leads to the expected high T 
behavior for the Polyakov loop when 1 < A < 4. 

It is clear that the allowed values of A become closer to 4 if a large number of higher-order 
derivatives of q vanish at the boundary. This should be expected since the case where all the 
derivatives of q are zero at the boundary corresponds to the exactly marginal deformation 
where A=4 and the resulting gauge theory is still a CFT where Uq and Cq are identically 
zero. It would be interesting to see how these consistency conditions change when multiple 
scalar fields and other supergravity modes are present such as in the supergravity dual of 
M = T theory [32] . 



The high T limit of Sq determines the value of I at the nontrivial UV fixed point, i.e, 

A hm /( T)=exp(|>)) (22) 

up to a T independent integration constant. When the conditions derived above are satisfied, 
£(T) approaches its asymptotic value in Eq. f l22l) from below in the class of gravity duals 
considered here. 

The nontrivial restrictions on the scaling dimension of scalar fields in gravity duals derived 
here may prove useful in practice to develop consistent deformations of CFT's that can 
be used as effective models for the strongly-coupled quark-gluon plasma {33 1 produced in 
ultrarelativistic heavy ion collisions. 

I am grateful to A. Dumitru for many helpful comments and M. Gyulassy, R. Pisarski, 
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